ON THE EXISTENCE OF SOLUTIONS TO THE OPERATOR RICCATI 
EQUATION AND THE TAN THEOREM 
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Abstract. Let A and C be self-adjoint operators such that the spectrum of A lies in a gap of 
the spectrum of C and let li > be the distance between the spectra of A and C. We prove that 
under these assumptions the sharp value of the constant c in the condition ||B|| < cd guaranteeing 
the existence of a (bounded) solution to the operator Riccati equation XA — CX + XBX = B* is 
equal to V2- We also prove an extension of the Davis-Kahan tan theorem and provide a sharp 
estimate for the norm of the solution to the Riccati equation. If C is bounded, we prove, in 
addition, that the solution X is a strict contraction if B satisfies the condition < d, and that 
this condition is sharp. 



1. Introduction 
In this paper we consider the operator Riccati equation 

XA-CX+XBX = B* (1.1) 
associated with the self-adjoint 2x2 block operator matrix 

"=(b- c) C-^* 

on the orthogonal sum IH' = ^H'a® of separable Hilbert spaces 9{a and 0-Cc- Here A is a 
bounded self-adjoint operator on the Hilbert space C possibly unbounded self-adjoint op- 
erator on the Hilbert space 9{c, and 5 is a bounded operator from 9{c to ^a- 

Solving the Riccati equation appears to be an adequate tool in the study of the invariant 
subspaces of the operator H that are the graphs of bounded operators from J-Ca to J-Cc, the so- 
called graph subspaces. It is well known that given a bounded solution X : 9{a to the 
Riccati equation (11.11) (with RanZ C Dom(C)), the graph 

g{x) = {x®Xx\ X e ^a] 

of the operator Z reduces the operator// (see, e.g., 01 Section 5]) and in the framework of this 
approach the following two problems naturally arise. The first problem is to study the spec- 
trum of the part of H associated with the reducing subspace g {X) (respectively its orthogonal 
complement Q {X)^), and the second one is to estimate the operator angle (see, e.g., flO| for 
discussion of this notion) between the subspaces Q (X) and 9{a (respectively g {X)^ and 0{c)- 
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Both of these problems can efficiently be solved if a bounded solution X to (11.11) is known: the 
operator H appears to be similar to the diagonal block operator matrix 

A+BX 
C-B*X' 

associated with the decomposition 9{ = 9{a® (in particular, o{H) = o{A + BX) U o(C - 
B*X*)), and the operator angle between the subspaces g {X) and 0-Ca has the representation 



However, Cbest is known to be in the interval 
then Cbest e 



= arctanVX*X. 

If the spectra of A and C overlap, the Riccati equation (11.11) may have no solution at all 
(cf., e.g., [3, Example 3.2]). At the same time the spectra separation requirement alone does 
not guarantee the existence of solutions either (see, e.g., [T, Lemma 3.11]). Under the spectra 
separation hypothesis 

dist(o(A),a(C)) >0, (1.3) 

a natural sufficient condition for the existence of solutions to the Riccati equation (11.11) requires 
a smallness assumption on the operator B of the form 

<Cbestdist(o(A),o(C)) (1.4) 

with a constant Cbest > independent of the distance between the spectra o(A) and o(C) of 
the operators A and C, respectively. The best possible constant Cbest in ( 11-41 ) is still unknown. 

TT^^ \/2 (see |l3|). If both A and C are bounded, 
cn, V2\ with Ct, = ^"~^^^'+^^ = 0.503288... (see O)- In O the best possible 

constant Cbest has been conjectured to be Some earlier results in this direction can be 

found in[2|, |17|,|19|, and|20|. 

In some particular cases the optimal solvability condition (11.41 ) can be relaxed provided that 
some additional assumptions upon mutual disposition of the spectra of A and C are posed. For 
instance, if the spectra of A and C are subordinated, e.g., 

supo(A) <info(C), (1.5) 

the Riccati equation (11.11 ) is known to have a strictly contractive solution for any bounded B 
(see, e.g., 1 1 1). To some extent abusing the terminology one may say that in this case the best 
possible constant in inequality (11.41) is infinite: No smallness assumptions on B are needed. 
In the limiting case of (11.51) . 

sup a(A) = inf o(C), 

the existence of contractive solutions has been established in fT\ under some additional assump- 
tions which have been dropped in lITTll . See also [18| where the spectra separation condition 
( 11.31) has also been somewhat relaxed and the existence of a bounded but not necessarily con- 
tractive solution has been established. 

Our first principal result concerns the case where the operator C has a finite spectral gap 
containing the spectrum of A. Recall that by a finite spectral gap of a self-adjoint operator T 
one understands an open finite interval on the real axis lying in the resolvent set of T such that 
both of its end points belong to the spectrum of T. 
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Theorem 1. Assume that the self-adjoint operator C has a finite spectral gap A containing the 
spectrum of the bounded self-adjoint operator A. 

( i) Suppose that 

||5|| < where J = dist(o(A),a(C)), (1.6) 

with I ■ I denoting Lebesgue measure on R. Then the spectrum of the block operator matrix 

H = ^ in the gap A is a (proper) closed subset of( the open set) A. The spectral subspace 

of the operator H associated with the interval A is the graph of a bounded solution X : 

to the Riccati equation (11.11) . Moreover, the operator X is the unique solution to the Riccati 

equation in the class of bounded operators with the properties 

o{A + BX)cA and o{C - B*X*) cR\ A. (1.7) 



( ii) If in addition the operator C is bounded and 

\\B\\ < ^/d{\A\-d), (1.8) 

then the solution X is a strict contraction and 

1 / 2||A5 + 5C|| \ , _ 
\\X\\ < tan-arctan ^ — < 1. (1.9) 

" 2 \d{\A\-d)-\\Bf J ^ ' 

As a corollary, under the assumption that the operator C has a finite spectral gap A containing 
the spectrum of A, we prove that c = Vl is best possible in condition (11.41) ensuring the existence 
of a bounded solution to the Riccati equation (11.11) (see Remark l33t while c — 1 is best possible 
in (11.41) to ensure that the solution is a contraction (see Remark 1531 . 

The proof of the part (i) of Theorem [l] will be given in Section |3] and that of the part (ii) in 
Section 13 

Our second principal result holds with no a priori assumption upon the mutual disposition 
of the spectra of A and C (in particular, the spectra of A and C may overlap). 

Theorem 2. Assume that the self-adjoint operator C has a spectral gap A (finite or infinite) and 
the self-adjoint operator A is bounded. Assume that the Riccati equation (11.11) has a bounded 
solution X and hence the graph subspace Q (X) reduces the block operator matrix H. Suppose 
that the spectrum of the part li\g(^x) of the operator H associated with the reducing subspace 
^ {X) is a (proper) closed subset of (the open set) A. Then the the following norm estimate 
holds: 

11-611 

Pll<^ w?Y/j 5 = dist(o(//|^(;^)),a(C)). (1.10) 



Equivalently, 



tan0||<H^ (1.11) 



where is the operator angle between the subspaces <^nd Q{X). 
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Estimate (11.1 II) extends the Davis-Kahan tan0 theorem Q, a resuk previously known only 
in the case where the spectra of C and subordinated, that is, the operator C is semi- 

bounded and the spectrum of the part H\g!^x) li^s in the infinite spectral gap of C. This general- 
ization extends the list of the celebrated Davis-Kahan sin0 and sin 20 theorems, proven in the 
case where the operator C has a gap of finite length ||2l|. 

The proof of Theorem |2l will be given in Section |4l 

Our main techniques are based on applications of the Virozub-Matsaev factorization theorem 
for analytic operator- valued functions [ |22| (in the spirit of the work JEI (cf. |15|)) and the 
Daletsky-Krein factorization formula [5|. Under the hypothesis of Theorem[T]we prove that 

• for X ^ o(C) the operator-valued Herglotz function M(k) = XI -A + B{C -Xl)^B* 
admits a factorization 

M{X) =W{X){Z-XI), (1.12) 

with W being an operator- valued function holomorphic on the resolvent set of the op- 
erator C and Z is a bounded operator with the spectrum in the spectral gap A of the 
operator C, 

• the Riccati equation (11.11 ) has a bounded solution of the form 

X = — - f dX{C-Xl)-^B*{Z-Xl)-\ (1.13) 
2711 Jr 

where F is an appropriate Jordan contour encircling the spectrum of the operator Z, 

• the spectral subspace of the 2x2 block operator matrix H (11.21 ) associated with the 
interval A is the graph of the operator X, 

• the spectrum of the operator H in the interval A coincides with that of the operator Z, 
thatis, o(//)nA = a(Z). 

In Section 121 we recall the concept of invariant graph subspaces for linear operators as well 
as their relation to the Riccati equation. Theorem 12 . 51 below presents a general result linking 
the factorization property (11.121) of the operator valued function M{X) with the existence of a 
spectral subspace for the 2x2 self-adjoint block operator matrix ( 11.21 ) admitting representation 
as the graph of the operator (11.131 ). In Section |3] under hypothesis (11.61 ) we prove factorization 
formula (11.121 ) and give bounds on the location of the spectrum of the operator Z (Theorem 
13.21 ). and finally prove the part (i) of Theorem[T] The proof of TheoremElis given in Sectional 
In Section |5] combining the results of Theorems [T] (i) and launder assumptions (11.61) and (11.81) 
respectively we provide norm estimates on the solution X of the Riccati equation and prove 
Theorem 1 (ii). We conclude the section by an example showing that condition (11.81) ensuring 
the strict contractivity of the solution X is sharp. 

Few words about notations used throughout the paper. Given a Hilbert space !}(_ by ly^^ we 
denote the identity operator on . If it does not lead to any confusion we will simply write / 
instead of more pedantic notation . The set of all bounded linear operators from the Hilbert 
space ?iC to a Hilbert space L will be denoted by S ( ?iC , ) • If = ^ the shorthand « ( ) will 
be used for this set. Let K and L be self-adjoint operators on a Hilbert space 3C ■ We say K < L 
(or, equivalently, L > ^) if there is a number Y> such that L — ^ > yl. The notation p(r) will 
be used for the resolvent set of a closed operator T. 

After completing this work we learned that a result similar to the part (ii) of Theorem[T]has 
been recently obtained within a different approach by A. V. Selin (private communication). 
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2. Invariant graph subspaces and block diagonalization 



In this section we collect some results related to the invariant graph subspaces of a linear 
operator as well as to the closely related problem of block diagonalization of block operator 
matrices. 

First, recall the definition a graph subspace. 

Definition 2.1. Let OChea closed subspace of a Hilbert space and X G S ( ^ , 3iC ^ ) • Denote by 
and Pj^± the orthogonal projections in fAC onto the subspace !Ki and orthogonal complement 
!]<C^, respectively. The set 

g{X) = {xe9l\ P^±x = XP^x] 

is called the graph subspace associated with the operator X. 

For notational setup we assume the following 

Hypothesis 2.2. Let Hq be a self-adjoint operator in a Hilbert space ^ and a reducing 

subspace o/Hq. Assume that with respect to the decomposition 



with A being a bounded self-adjoint operator in ^a, C a possibly unbounded self-adjoint op- 
erator in ^c> cind Dom(//o) = © Dom(C). Assume, in addition, that with respect to the 
decomposition (12.11) the self-adjoint operator H reads as 



where B is a bounded operator from to ^a- 

Under Hypothesis 12.21 a bounded operator X from ^H^a to is said to be a solution to the 
Riccati equation (11.11) if Ran (X) C Dom(C) and (11.11) holds as an operator equality. 

The existence of a bounded solution to the Riccati equation (11.11) is equivalent to a possibility 
of the block diagonalization of the operator matrix H with respect to the decomposition ^ — 
g (X) © g (X)^. The precise statement is as follows (see Lemma 5.3 and Theorem 5.5 in [3|; 

also cf. 0, la, EDI, EOI). 

Theorem 2.3. Assume HyDothesis \2.2\ Then a bounded operator X from "Ha to Hq is a solution 
to the Riccati equation (11.11 ) iff the graph g (X) ofX reduces the 2x2 block operator matrix H. 
Moreover, ifX G (B {:Hat^c) ^ solution to (11.11) then: 

(i) The operator V^HV with 



the operator Hq reads as the block diagonal operator matrix 

//o = diag(A,C) 



(2.1) 




(2.2) 




is block diagonal with respect to decomposition (12.11) . Furthermore, 
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where Z = A + BX with Dom(Z) =^KAandZ = C- B*X* with Dom(Z) = Dom(C). 
(ii) The operator 

A={I + X*Xy/^Z{I + X*X)-^^^ (2.3) 
and possibly unbounded operator 

A={I + XX*)^/^Z{I + XX*)-^/^ (2.4) 

with Dom(A) = (/ + XX*)^/^(Dom(C)) are self-adjoint operators in J-[a <^nd 9{c, respectively. 

Theorem 12.31 yields the following uniqueness result as a corollary. 

Corollary 2.4. Assume Hypothesis \2.2\ Suppose that E and E are disjoint Borel subsets o/M 
such that dist(E,E) > 0. Let X = X (A,5,E,E) be the set of all bounded operators X from J-[a 
to 9-[c with the properties 

a(A + 5Z) cE, (2.5) 

o(C-5*X*)cE, Dom(C-5*X*) =Dom(C). (2.6) 

Then ifX, 7 G X satisfy the Riccati equation (11.11 ). then X = Y. 

Proof. Suppose X and Y are two bounded solutions to (11.11 ) both satisfying ( 12.51) and ( 12.61 ). Then 
by Theorem 12.31 the graphs of X and Y both coincide with the spectral subspace of the 2x2 
operator matrix ( 12.21) associated with the set E, and hence, X = Y . □ 

Under Hypothesis 12 . 21 introduce the operator-valued Herglotz function 

M{1)=XI-A + B{C-XI)-'^B\ Xep(C). (2.7) 

By definition the spectrum a(A/) of the function M is the set of all ?i G C such that either the 
operator M(?i) is not invertible or the inverse [M(A,)]^^ is an unbounded operator. 

It is well known (see, e.g., 1 17 1) that the resolvent of the operator H (12.21) can be represented 
as the following 2x2 operator matrix 

'0 

^0 (c-x/)-i 

-( _(C-L)-i5* )^(^)"'(^ -5(C-X/)-i), (2.8) 
Xep(//), 

where M is the Herglotz function given by (12.71 ). Representation ( 12.81) shows that for X E p(C) 
the operator H — XI has a bounded inverse iff M{X) does, which means that 

o(//)np(C) = o(M)np(C). (2.9) 

We will also need the following general result (cf. [2, Theorem 2.2], ifT^ Proposition 2.4 and 
Theorem 2.5], and [|15t Theorems 4.4 and 5.1]). 

Theorem 2.5. Assume Hyvothesis \2.2\ and suppose that the Herglotz function ( 12.71 ) admits the 
factorization 

M{X)=W{X){Z-XI), Xe^, (2.10) 
where Z is a bounded operator in such that o(Z) fl o(C) = 0, Q is a domain in p(C) such 
that (5 (Z) C and W is a holomorphic 'B (^N^a) -valued function on Q, such that for any X G CJ(Z) 



{H-XI) 
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the operator W{X) has a bounded inverse. Then o{Z) is an isolated part of the spectrum of the 
operator H and the spectral subspace RanE/f (o(Z)) ofH associated with the set o{Z) is the 
graph of the bounded operator X from to given by 

X = -^ I dXiC~Xiy^B\Z-Xl)-\ (2.11) 
zra Jt 

Here, T is an arbitrary Jordan contour in p (Z) fl p (C) ( maybe consisting of several simple 
Jordan contours) encircling o(Z) in the clockwise direction and having winding number with 
respect to the spectrum ofC. 

Moreover, the operator Z can be written in terms of the operator X given by ( I2.11I) as 

Z = A + BX (2.12) 

and the factor W{X) admits analytic continuation to the whole resolvent set of the operator C 
by the following formula 

W{X)=I-B{C-X)-^X, lep{C). (2.13) 

Proof. By hypothesis the function W{X) is holomorphic on an open set Q. C p(C) containing 
the closed subset o(Z) and the operator W{X) has a bounded inverse for any X E o(Z). Hence 
there is an open neighborhood Q. of o(Z) in where the operator W{X) is boundedly invertible, 
i.e. W{X)-^ E 'B {^a) for any lEQ..By (ITTIl 

M{X)-^ = {Z-Xl)-^W{X)-\ Xea\o(Z). (2.14) 

Taking into account (12.91) one infers that the spectrum ofH'mQ. coincides with that of M and, 
thus, with that of Z, that is, 

o(//)na = o(z), (2.15) 

Since Q. is an open set, o(Z) is a closed set, and o(Z) C Q., one also concludes that a(Z) is 
isolated from the remaining part of the spectrum of H, i.e. 

dist(a(Z),a(//)\a(Z)) >0. 

Using representation ( 12. 8I ). for the spectral projection Eh (o(Z)) of the operator H associated 
with the set o(Z) the Riesz integration yields 

E^(a(Z)) =E^(o(Z)) = dX{H-Xl)-' 

r 

= -Si/''K-(c-W)-'r)^W' -s(c-W)->), 

where T stands for an arbitrary Jordan contour (possibly consisting of several simple Jordan 
contours) in Q. encircling the spectrum of Z in the clockwise direction and having winding 
number with respect to the spectrum of C. Hence, 




(2.16) 
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where 



E 
F 



/ dX{C -XI)-^B*M{X)-^B{C -Xl)-\ 
2ra Jr 



(2.17) 



and 



v-l 



v-l 



G=-L / dX{C-Xl)-^B*M{X) 
2ni Jr 

= ^ I dX{C-Xl)-^B*{Z-Xl)-^W{X)-\ 

using factorization formula (I2.10I) . Since both the operator-valued functions (C — X)^^ and 
W{X)^^ are holomorphic in Q., applying the Daletsky-Krein formula (see |5 , Lemma 1.2.1]) we 
get 

/ dX{C-Xl)-^B*{Z-Xl) 
2ni Jr 

/ dX{Z-Xl)-^W{X) 
2711 Jr 

Hence, combining (12.141) and (12.171) proves the representation 

G = XE, 

where X (the first factor on the r.h.s. part of (12.181) ') is given by 

X = --^ / dX{C-Xl)-^B\Z-Xl) 
Ini Jr 

In an analogous way one also proves that 

F=XEX*. 

Clearly, for XeT we have Ran(C — Xl)^^ C Dom(C) and hence 

Ran(X) C Dom(C), 

which immediately follows from (12.201) . Multiplying both sides of (12.201) by B from the left 



(2.18) 
(2.19) 
(2.20) 
(2.21) 



yields 



5Z = --^ / dXB{C-Xl)-^B\Z-Xl)-\ (2.22) 
2711 Jr 

Meanwhile, 

B{C-XI)^B* =A-XI+M{X) =A-XI + W{X){Z-XI), XeT, 
and, hence, using (12.221) 

BX = -^ [ dX\W{X)+A{Z-Xiy^-X{Z-Xl)-^]. (2.23) 
2711 Jr 

The function W{X) is holomorphic in the domain bounded by the contour F and, thus, the first 
term in the integrand on the r.h.s. of (12.231) gives no contribution. Since F encircles the spectrum 
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of Z, the integration of the remaining two terms in (12.231) can be performed explicitly using the 
operator version of the residue theorem, which yields BX = —A + Z and hence 

Z = A + BX (2.24) 

proving (12.121 ). 

Since the spectra of the operators C and Z are disjoint and Z is a bounded operator, it is 
straightforward to show (see, e.g., |4| or [21 1) that the operator X given by ( 12.201) is the unique 
solution to the operator Sylvester equation 

XZ-CX = B\ 

which by (12.241) proves that X solves the Riccati equation dl.ll) . 

Now applying Theorem 12.31 and using (12.161) . ( 12.191 ). and (12.211) we arrive at the series of 
equalities 



^:^J^dXPAV-\H-Xl)-WPl 
^PAV-'EH{oiZ))VPl 

^E{I + X*X), (2.25) 



where Z = C-B*X* with Dom(Z) = Dom(C), 

V 



I -X' 
X I 

and Pa is the canonical projection from ^M' to iM'a (i.e., PaHa ® fc) = /a ior /a E ^H^a and 
fc £ ^c)- Combining ( 12.191 ). and (12.251) one concludes that the spectral projection E/f (a(Z)) 
admits the representation 

^H{o{Z)) - [^x{i + x*x)-' xii + x*x)-'x* J ■ ^^-^^^ 

Note that the contour F in ( 12.201 ) can be replaced by an arbitrary Jordan contour in p (Z) fl p (C) 
(possibly consisting of several simple Jordan contours) encircling the set o(Z) in the clockwise 
direction and having winding number with respect to o(C). Then observing that the r.h.s. of 
(12.261) is nothing but the orthogonal projection in i?/" = i?/^ © fM'c onto the graph of the operator 
X proves that Ran E//(a(Z)) = g{X). 

We conclude with the proof of the representation (12.131) . First, by (12.241) we notice that 

[I-B(C-iy^X]{Z-X) = [I-B{C-X)-^X]{A+BX-'k) 

= A+ BX - X- B{C - XY\XA + XBX - XX) , 
p(C). 



10 V. KOSTRYKIN, K. A. MAKAROV, AND A. K. MOTOVILOV 

Since X solves the Riccati equation (11.11) one infers that XA + XBX = B* -\- CX which implies 

[I-B{C-Xy ^X] {Z~X)=M{X)+BX-B{C-Xy\C- X)X 

= M(X), XGp(C). (2.27) 

Hence combining (12.101 ) and (12.271 ) yields 

w{X) =M{X){z-XY^ =i-B{c-XY^x, iea\o{z). 

Then the analytic continuation completes the proof. □ 

Remark 2.6. By analytic continuation of both parts of (12.101 ) one concludes that the factorization 
formula (12.101 ) holds for any X e p(C) with W{X) given by (12.131 ). Moreover, the representation 

[W{X)]-^ = (Z - XI) [M{X)]-\ Xea\ a{M) = Q. \ a(Z), 

and the fact that o(M) C o{H) imply that admits analytic continuation as a 'B{^a)- 

valued holomorphic function to the domain p{H) U o(Z). 

3. An existence result. Proof of Theorem 1 (i) 

As we have already mentioned in Introduction our main technical tool in proving the solv- 
ability of the Riccati equation is the Virozub-Matsaev factorization theorem ll^ (also see [ 16 1). 
For convenience of the reader we reproduce the corresponding statement following Propositions 
1.1 and 1.2in [18|. 

Theorem 3.1. Let ^ be a Hilbert space and F{X) a holomorphic 'B {nc) -valued function on a 
simply connected domain Q C C. Assume that Q. includes an interval [a,b] C M such that 

F(a)<0, F{b)>0, and ■^F{X)>0 for all Xe[a,b]. 

dX 

Then there exist a domain Q.GQ. containing [a, b] and a unique bounded operator Z on with 
o(Z) C {a,b) such that F{X) admits the factorization 

F{X) = G{X){Z-XI), Xeh, 

where G{X) is a holomorphic operator-valued function on Q. whose values are bounded and 
boundedly invertible operators in Ki, that is, 

G{X)e'B{'K) and [G{X)Y^ e 'B{'K) Xeh. 

Based on Theorem 13. II we obtain the following factorization result, the cornerstone for our 
further considerations. 

Theorem 3.2. Assume Hypothesis \2.2\ Assume, in addition, that C has a finite spectral gap 
A — (a, (3), a < P, the spectrum of A lies in A, i.e., o(A) C A, and 

\\B\\ < ,/d\Al 

where 

J = dist(o(A),a(C)). 

Then there is a unique operator Z G'B (^^a) v^^ith o(Z) C A such that the operator-valued func- 
tion M(X) given by (12.71) admits the factorization 

M{X)=W{X){Z-XI), Xep{C), (3.1) 
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with a holomorphic 'B {9{a) -valued function W on p(C). Moreover, for any 

Xe (C\a(M))uA 
the operator W (k) has a bounded inverse and 

a(Z) =o(//)nAc [infa(A) - 5_, sup a(A) +5+], (3.2) 

where 

-arctanTT — . ' , , ) < info(A) - a, (3.3) 
2 p — info(A) / 

/I 2 Il5 II \ 
5+ = ||5||tan I -arctan I <P-supo(A). (3.4) 



2 supo(A)-a 
Proof. By the spectral theorem 

B{C-II) 'B*= I BEc{dM)B*j^. ^ep(C), 

K\A 

where Ec(ju) stands for the spectral family of the self- adjoint operator C. Hence 



(3.5) 

For Xe A the integral in ( 13.51) is a non-negative operator. Therefore, the derivative of M(X) is a 
strictly positive operator: 

-^M(?i) > / > 0, leA. 
dk 

Next we estimate the quadratic form of M(?i). Let / G 0-Ca, \\f\\ = 1- Then 

(M(X)/, /) = X - (A/, /) + ( (C - X) - 57) 

a 

= X-(A/,/) + j -l^{Ec{dij)B*fXf) (3.6) 

— 00 

-|-oo 

+ / ^(Ec W57,57), X G p(C). 

Since for ?i G A the integral in the second line of ( 13.61) is non-positive and the one in the third 
line is non-negative, one obtains the two-sided estimate 

ll5lP\ / Il5ll^\ 

;i- sup a(A) - /<M(;i) < X - inf a(A) - /, X G A. 



k-aj ~ V 
Now, a simple calculation shows that 

M{k)<0 for Xg (a,info(A)-6_) (3.7) 

and 

M{k)>0 for Xg (sup o(A) + 6+, p), (3.8) 
where 5- and 5+ are given by (13.31) and (13.41) . respectively. 



12 



V. KOSTRYKIN, K. A. MAKAROV, AND A. K. MOTOVILOV 



Thus, the function F(?i) = M{X) satisfies assumptions of Theorem B . 1 I for any a G (a, inf a(A) - 
6_) and any b E (supa(A) +5+,|3), proving the existence of the unique bounded operator 
Z e 'B {^Ha) such that O.ll) and (13.21) hold taking into account ( 12.91) . It follows from Theorem 
13. II that the factor W{X) in (12.101) has a bounded inverse in a complex neighborhood t/ C C of 
the interval [inf a(A) — 5-, sup a(A) +5+]. Moreover, the operator W{X) has a bounded inverse 
for any ?iG (C\o(M))uAby Remark I231 The proof is complete. □ 

We are ready to prove the part (i) of Theorem 1 . 

Proof of Theorem\I\(i). By Theorem 13.21 the Herglotz function ( 12.71) admits the factorization 
(12.101 ) with W(k) and Z satisfying hypothesis of Theorem 12.51 Therefore, the Riccati equation 
(11.11) has a bounded solution X given by (12.1 II) . Theorem l2.5l also shows that the graph g (X) of 
the operator X coincides with the spectral subspace Ran Eh (<5(Z)) for the block operator matrix 
H given by (12.21 ). By ( 13.21 ) the subspace Ran (<3(Z)) coincides with Ran E^ (A) , proving that 

RanE/f(A) = ^(X). (3.9) 
To prove thatX possesses the properties (11.71) we proceed as follows. Let 

^=(^X ~f*)' S={I + X*Xy/^, and S={I + XX*y/^. (3.10) 
From Theorem l2.3l it follows that 

where A and A are the self-adjoint operators defined by (12.31) and (12.41) . respectively; Ea(A) and 
E^(A) denote the spectral projections for A and A associated with the interval A. Since A is 
similar to Z (A = SZS^^) and by Theorem 13. 21 the inclusion o(Z) C A holds, one concludes that 
o(A) C A. Hence, Ea(A) = / and then (13.111) implies the equality 

/ \ , „ / \ 



E.(A)=y( oj^ +^0 5-iE;^(A)5 )^ ■ ^''''^ 



The first summand on the r.h.s. of (13.121 ) 



/ \ I _ f {i + x*x)-^ {i + x*x)-^x* 



^yo )^ \x{i+x*x)-^ x{i+x*x)-^x* 

coincides with the orthogonal projection onto the graph subspace g (X) and, hence, by (13.91 ) it 
equals Eh (A) . Thus, the second summand on the r.h.s. of (13.121) vanishes. 







which means that E^(A) = and, therefore, o(A) fl A = 0. By (12.41) it is now straightforward 
to see that X satisfies (11.71) . It remains to conclude that by CoroUarv 12.41 the operator X is 
the unique bounded solution to (11.11) satisfying (11.71) which completes the proof of Theorem 1 
(i). □ 
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Remark 3.3. Obviously, under hypothesis of Theorem [T] the inequality | A| > 2d holds, with the 
equality sign occurring only if the spectrum of the operator A is a one point set. Hence, the 
condition ||5|| < ^/2d, which is stronger than (11.61) . implies the existence of a bounded solution 
to the Riccati equation dl.ll) . This means that in the case where C has a finite spectral gap A and 
o{A) C A the best possible constant Cbest in condition (11.41) ensuring the solvability of the Riccati 
equation satisfies the inequality Cbest > V^. On the other hand in |T, Lemma 3.11 and Remark 
3.12] it is shown that Cbest < V2. Thus,c = a/2 is best possible in inequalitv (II .41) which ensures 
the solvability of the Riccati equation under the additional hypothesis that dist(o(A) , o(C) ) > 
and that the spectrum of A lies in a spectral gap of the operator C. 

4. The tan© Theorem 

We start out by recalling a concept of the operator angle between two subspaces in a Hilbert 
space going back to the works by Friedrichs [8|, M. Krein, Kransnoselsky, and Milman (13], 
iri4J . Halmos ||9J, and Davis and Kahan [TJ. A comprehensive discussion of this concept can be 
found, e.g., in ifTOl . 

Given a closed subspace Q, of the Hilbert space = 9{a® introduce the operator angle 
between the subspaces © {0} and Q. by 

= arcsin^/^^-PAQPl, (4.1) 

where is the canonical projection from ^ onto and Q the orthogonal projection in ^ 
onto (Z- If the subspace Q, is the graph (X) of a bounded operator X from ^a to iHc, then (see 
ifTOl : cf. and 0) 

tm® = Vx*X (4.2) 

and 

||sin0|| = ||Q-P||, (4.3) 

where P = Pa* Pa denotes the orthogonal projection in ^ onto the subspace ^a © {0}. 

Note that the common definition of the operator angle (see, e.g., fTOl) slightly differs from 
(14.11) . Usually, the operator angle is defined as the restriction of the operator (14.11) onto the 
maximal subspace of ^a where it has a trivial kernel. Clearly, the difference in these two 
definitions does not effect the value of the norm || tan0||. 

Now we are ready to prove the second principal result of the paper, a generalization of the 
Davis-Kahan tan0 Theorem Q. 

Proof of Theorem^ By hypothesis the Riccati equation dl.ll) has a bounded solution X. Then, 
by Theorem l2.3l the operator Z = A + BX is similar to the bounded self-adjoint operator A given 
by (12.31 ) and hence 

o(Z)=a(A)cM. (4.4) 
The Riccati equation (11.11 ) can be rewritten in the form 

X{Z-yI)-{C-yI)X = B\ (4.5) 

where 

Y = ^(supa(Z)+infa(Z)) = ^ (supo(A) +infa(A)). (4.6) 
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By hypothesis dist(a(//|^(jc)),a(C)) = 5 > and then Theorem 12 . 31 implies that 

<7(^l^w)=o(Z)=o(A), 
which by hypothesis proves the inclusion 

o(C) C (-°o,info(A) -6] U [supo(A) +5,°o). (4.7) 
Hence, combining (14.41) . (14.61) . and (14.71) proves that y G p(C) and 

"(^-«"« = F^^- <^-^' 

Multiplying both sides of (14.51) by (C — y/)^^ from the left one gets the representation 

X = {C-yiy\x{Z-yI)-B*). (4.9) 

Using Theorem 12. 31 (11) one obtains the estimate 

||X(Z-y/)|| =\\X{I+X*X)-^/^{A-yI){I + X*X)^/^\\ 

< \\X{I+X*X)-^/^\\\\A-yI\\{l + \\Xfy/^. (4.10) 

Clearly, by the spectral theorem 

\\X\\ 



||x(/+x*x)-i/2|| = ^||x*x(/+x*x)-i|| = -^ll^l|2)i/2- 

Hence (14.101) implies the estimate 

||X(Z-y/)||< 11X11 l|(A-y/)l|, 
which together with ( 14.91) proves the norm inequality 

P||<||(C-y/)-i||(||A-y/||||X|| + ||5||). (4.11) 

Solving inequality (14.1 II) with respect to ||X|| and taking into account (14.81) proves the norm 
estimate (11.101 ). Finally, since || = by the definition of the operator angle ( 14.21) one 

gets 

||tan0|| = 

Hence, dl.lOl) is equivalent to (11.111) . □ 

Remark 4.1. It is natural to ask whether estimate dl.l II ) remains to hold if one replaces the 
distance 5 = dist(a(//|^(;c)),a(C)) by 

5 = dist {o{H\ g (;^)x , o(A) ) = dist (a(C - B*X*), o(A)) . 

The answer is "No": Example 6.1 in |7J shows that the inequality 5|| tan0|| < fails to hold 
in general. 
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5. Norm Estimates of Solutions. Proof of Theorem 1 (ii) 

The existence result of Theorem 1 (i) by itself gives no clue for estimating the norm of the 
corresponding solution X to the Riccati equation. To the contrary, Theorem 2 provides such 
an estimate whenever some additional information on the spectrum location of the "perturbed" 

operator matrix H = is available. In turn, the bounds on the spectrum of the part 

H\g(^x) of the operator matrix H associated with the reducing subspace g (X) (needed to satisfy 
the hypotheses of Theorem |2l) can be obtained by combining the results of Theorems 12.51 and 
13.21 As a result of performing this program one gets an a priori estimate on the norm of the 
solution X. 

Theorem 5.1. Assume hypothesis ofTheorem\^(i) with A = (a, (3), a < (3. LetX be the unique 
solution to the Riccati equation (11.11) referred to in Theorem\I\ Then 

||X||<H, (5.1) 
5 

where 

5 = min{info(A)-a-5_, p - supo(A) - 5+} > (5.2) 
with d± given by (13.31 ) and (13.41 ). 

Proof. Under hypothesis of Theorem [T](i) with A = (a, |3), a < (3 one can apply Theorem 13.21 
and using the same strategy of proof as that of Theorem[T](i) one concludes that 

Z=A+BX, 

where Z is the unique operator with o(Z) C A referred to in Theorem 13.21 and X is the unique 
solution to the Riccati equation (11.11) referred to in Theorem [l](i). Hence, by Theorem 13. 21 



o{A + BX) c [infa(A)-5_,supo(A)+5+]. (5.3) 
By hypothesis (of Theorem [T](i)) 

o(C)c (-oo,a]u[p,oo), (5.4) 

which together with (15.31) yields the inequality 

dist(o(A + 5X),o(C)) >5. 

By Theorem l2.5l one observes that o{H) \g(^x) = ^(2) which proves (15.11 ) using Theorem|21 □ 

The proof of Theorem [T](ii) needs complementary considerations. Our reasoning is based on 
the celebrated Davis-Kahan tan20-Theorem [7|. For convenience of the reader we reproduce 
the corresponding result (cf. Theorem 2.4 (iii) and Remark 2.8 in |11| and Corollary 6.4 in 
flOl). 

Theorem 5.2. Assume HyDothesis \2.2\ Suppose that the operator C is bounded and supa(A) < 
info(C). Then the open interval (supG(A),infa(C)) is a spectral gap of the operator H and 
the spectral subspace E// ((— oo, supo(A)]) is the graph of a contractive operator X from 9{a to 
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:Hc- Moreover, the operator X is the unique contractive solution to the Riccati equation (11.11) 
and its norm satisfies the estimate 

n n 2||5||\ 

\\X\\ < tan -arctan < 1, 

where 

J = dist(o(A),o(C)). 

Now we are prepared to prove Theorem 1 (ii). 

Proof of Theorem\I\(ii). Let X be the solution to the Riccati equation (11.11 ) referred to in The- 
orem [l](i) and thus the spectral subspace of the operator H associated with the interval A is 
the graph <^{X) of the operator X. Then ^ (X) is also the spectral subspace of the operator 
{H — yl)^ associated with the interval [0, |Ap/4) where y is the center of the interval A, that is, 

^(X) =Ran(E^(A)) =Ran(E(^_Y/)2[0,|A|V4)). (5.5) 

By inspection one obtains that with respect to the decomposition fH' = 9{a®^a the non-negative 
operator {H — yl)^ reads 

iH-yi)'=(i fy (5.6) 



,5* 

whereA = {A-yI)^ + BB*,B=AB + BC, and C = {C -yl)^ +B*B. 

The hypothesis that the spectrum of C lies in M \ A implies the operator inequality 

- |Ap , lAp 
C> ^-^I + BB* > ^I. 
~ 4 ~ 4 

The hypothesis o(A) C A yields 

o I A 



0<A< ( ^-J 1 I + BB* < 



2 



d] +\\B 



|2 



(5.7) 



taking into account that dist(o(A) , o(C)) = d. Hence, under hypothesis dl.SI) one concludes that 

dist(a(A),a(C)) > d{\A\-d) - \\Bf > (5.8) 

and that the spectra o(A) and o(C) of the entries A and C are subordinated, that is, sup a(A) < 
inf o(C). By Theorem l5.2l (cf. Theorem 2.1 in |[T]|) one infers that the interval (supo(A),infa(C)) 
lies in the resolvent set of the operator [H — yl)^. In particular, the following inclusion holds 

((|A|/2-J)2+||5||2,|A|V4) cp((//-Y/)2). (5.9) 

Therefore, the spectral subspaces of the operator [H —ylY associated with the intervals [0, |Ap/4) 
and [0, |(A|/2 — (i)2+ respectively, coincide, that is, 

£^(X)=Ran(E(^_,,)2([0,|A|2/4))) 

= Ran(E(^_,,p([0,(|A|/2-J)2-M|5f])). 

From (15.61) one concludes that the operator matrix (H — yl)^ is an off-diagonal perturbation of 
the matrix diagjAjC} diagonal with respect to the decomposition 0-[ = 0-Ca ® ^c- Applying 
again Theorem 15 .21 proves that the spectral subspace of the operator {H — yl)^ associated with 
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the interval [0, |(A|/2 — + is the graph of a contraction X satisfying the norm-estimate 
(I1.9I) . By (15.51 ) and (15.101) this subspace coincides with g {X) and, therefore, X = X and hence 
(11.91) holds for X, completing the proof. □ 

Remark 5.3. Condition (11.81) ensuring the strict contractivity of the solution X is sharp. This 
can be seen as follows. Let 9{a = C, 0-Cc = C^, A = 0, 

and 

By inspection one proves that the 2x1 matrix 




(5.11) 



solves the Riccati equation 

XA-CX+XBX = B*. 

Moreover, 

A+BX = 



and X possesses the properties (11.71) . Clearly, ||5|| = b and \\X 




For b G [d, \/2d) hypothesis (11.61) is satisfied with A = (—J, d), condition (11.81) fails to hold, 
and > 1, that is, estimate ( 11.81 ) is sharp. 
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